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We discuss Casimir effect of a massless, minimally coupled scalar field in a 6D warped flux 
compactification model and its implications for the hierarchy and cosmological constant problems, 
which are longstanding puzzles in phenomenology and cosmology. Due to the scale invariance of the 
background theory, the 4D effective theory contains a volume modulus. To stabilize the modulus, we 
consider one-loop corrections to the effective potential by the Casimir effect. The one- loop effective 
potential for the volume modulus has a form which is very similar to Coleman- Weinberg potential. 
We evaluate coefficients appearing in the effective potential by employing zeta function regularization 

t , and heat kernel analyses. The volume modulus is stabilized for smaller degrees of warping, below a 

critical value, which depends on deficit angle of the reference brane. After stabilizing the modulus, 
it is possible to obtain observed values of the mass ratio between the fundamental energy scales and 

s ! ' a tiny effective cosmological constant (though its sign is negative). The degree of warping should 

be tuned to be close to the critical value, not as severely as the original fine-tuning. 



o 
o 
o 



- 1—1 

X 



PACS numbers: 04.50.+h; 98.80.Cq 

Keywords: Extra dimensions, Quantum field theory 



o 

0" 1 I. INTRODUCTION 

U 

There are several longstanding problems in phenomenology and cosmology. One of them is why gravity is so weak 
in comparison with the electroweak interactions (in other words why Planck scale Mpi ~ 10 19 GeV is much larger than 
that of electroweak interaction Mew ~ 10 3 GeV) and is known as the hierarchy problem. Another important problem 
^ ' is why the energy density of dark energy which dominates the present Universe (assuming that its origin is vacuum 
CO energy of quantum fields) pvac ~ (10 3 eV) 4 is much smaller than that of gravitational scale Mpj, which is naturally 
expected from the standard model and is known as the cosmological constant problem. 

In this article, we focus on the Casimir effect in a 6D braneworld model whose extra-dimensions are compactified 
by magnetic flux 6D flux compactification models have attracted much attensions because of the fascinating 
feature that they may help resolve the above puzzles, namely the hierarchy problem [2] 1 and the cosmological 
constant problem [4J . It is also expected that in these models effects of extra-dimensions may be detected in future 
experiments on gravity because they may become important at sub-milimeter scales , where Newton's law has not 
been strictly confirmed. 

The basic motivation to consider the Casimir effect in such a 6D brane model is as follows. An important point 
is that Mpi and Mew have satisfy Mpi ~ 10 16 Mew and Py* c ~ 10~ 16 Mew- Thus in these two hierarchies there is 
substantial similarity and it could be naturally expected that if there is a theory which gives a common factor of 10 16 , 
both hierarchy problems may be solved at the same time. The Casimir effect in 6D spacetime may be able to give such 
a common factor. 2 Let's imagine a spacetime with n extra dimensions (later we set n = 2) whose size is assumed to 
be stabilized at a characteristic scale a. Then, the dimensionaly reduced Planck mass is given by Mpj = a n M%^™ 0], 
where we assume M± +n ~ Mew> Casimir energy density induced on the brane by fields living in the internal space 
is roughly given by peas ~ a -4 - Thus, we get 



M\ 8 /" 
, 4+n X 



(1.1) 



Especially, for the case of n = 2, we obtain pcas/M| ~ (Mb / Mpi) 4 . If the Casimir energy density peas plays the role 



* Email: masato" at" theorie.physik.uni-mucnchcn.de 

1 For extensions of the model given in with more successful localizations of the standard model particles on intersections of D-branes, 
see e.g., ||. 

2 After the proposal of Randall-Sundrum braneworld model Q, the Casimir effect in 5D braneworld has also been discussed in terms of 
phenomenology, especially toward a resolution of the hierarchy problem, see e.g., [f|. 
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of the dark energy density pvac, we can get the desired ratio (see also [7(). The problem is whether the ratio is really 
obtained from the setup of the 6D braneworld. Bearing the above considerations in mind, in this article, we focus on 
the Casimir effect in a specific model of 6D braneworld with a warped flux compactification. 3 

We analyze the Casimir effect induced by a massless minimally coupled scalar field in a warped flux compactification 
model based on a 6D (Salam-Sezgin) supergravity employing zeta function regularization techniques. In the 4D 

effective theory of the model, a volume modulus arises due to the scale invariance of the 6D theory and is stabilized by 
one-loop quantum corrections to the effective potential induced by the bulk scalar field. We then discuss implications 
for the hierarchy and cosmological constant problems. In Ref. [ll|, [T§] the Casimir effect in a 4D toy model, which 
has a very similar spacetime structure to this 6D model, has been discussed. As a consequence, it has been suggested 
that for larger degrees of warping and smaller degree of dilaton coupling, the Casimir effect may give a resolution to 
the hierarchy and cosmological constant problems. So, it is desirable to discuss the Casimir effect in the original 6D 
model. In attempting to evaluate the Casimir effect in the original 6D model, however, there is a technical problem 
arising from a lack of mathematical formulation of conical heat kernel coefficients (in 6D), which are relevant to 
contributions of the boundary 3-branes to the effective potential. 

To overcome the above technical problem, we instead focus on a special relation between the relevant heat kernel 
coefficient and the analytically continued zeta function [131 ]. The mass spectrum which includes brane contributions 
is now available in a conformally related (unwarped) spacetime and is very similar to the case of the 4D model. Thus, 
instead of trying to evaluate the relevant heat kernel coefficient, we shall use the analytically continued zeta function 
by employing the mode sum technique, which has been developed in the recent work on the 4D toy model, Ref. 

The article is organized as follows. In Sec. II, we briefly introduce the background model for a 6D warped flux 
compactification. In Sec, III, we discuss the quantum mechanical perturbations of a massless minimally coupled scalar 
field. In Sec. IV, we derive mathematical relations in order to analyze the stability of the volume modulus and the 
Casimir effect. In Sec. V, we will discuss volume stabilization and its implications for the Casimir effect on the 
hierarchy and cosmological constant problems. In Sec VI, we shall close this article after giving a brief summary and 
discussions on possible extensions of the present work. 



II. A WARPED CODIMENSION TWO BRANE MODEL WITH FLUX COMPACTIFICATION 



A. Solution 

We consider a 6D Einstein-Maxwell-dilaton theory with a non- vanishing scalar potential p], [l(| as 

S 6 = Ml J d 6 x^ Qi? - ^d A <pd A ip - \e^F AB F AB - 2g 2 e^ , (2.1) 

where <p is a dilaton field, Fab represents a U(l) gauge field strength and g is the UJX) gauge coupling constant. This 
theory corresponds to the bosonic part of the Salam-Sezgin, 6D supergravity [9|, [l0( . Hereafter we set Mi = 1 for 
simplicity and if needed, we put it back in explicitly. 

This theory contains a series of solutions of warped flux compactification [l[ ; 

do 2 

ds 2 = h{p)dd 2 + -f- + {2p)i 1 ^dx ti dx v , 

Hp) 

Hp) - £ (pI - p 2 ) (p'-p 2 -), 

p(p) = -ln(2p), 

p 9P+P- r9 9 x 

6p ~ p3~ ' ' 

where two 3-branes are located at p — p±. For later convenience, we define a new parameter 

a = 9 —, (2.3) 
P+ 



In the above discussion here, we implicitly assume that the tree level vacuum energy density is effectively cancelled by some self-tuing 
mechanism .4), though there are some criticisms for this mechanism We just discuss the impact of the Casimir effect (in a 6D model) 
on the cosmological constant problem and implicitly assume that a self-tuning mechanism exists. 



3 



which controls shape (i.e., warping) of the internal space. 

To see the spacetime structure in details, it is useful to introduce a new coordinate as 



P 2 +-P 2 \ 1/2 _5 2 (l-a 2 ) 

n2 _ „2 



V=—^ -0, (2-4) 



where n — 1 — 

Z7T 



p 2 — p_J 2k 



+ z 2 27T~T^2 m +V^^" 



g 2 (l + z 2 ) 2 g 2 (l + a 2 z 2 ) 2 



(2.5) 



S± = - d 4 xVha± , (2.6) 



ds 2 = (2p) 
The braneworld action is given by 



respectively, where a± denotes the brane tensions, which are related to to the conical deficit angles by 

a± = M£S± . (2.7) 

h a b is the brane induced metric. In order to share the same angle period around both the branes, the deficit angles 
are related to a as 

—a 2 - 2.8 

Ztt — 0_ 

Eq. (|2.8|) implies that once the brane tensions, a± are specified, then the bulk shape a is also fixed. We now regard a 
and k as free parameters, instead of <r±, along with g. So we shall use (+)-brane as a reference brane. The remaining 
modulus is the absolute size of the bulk, p + . Note, however, that there is also magnetic flux constraint given by 

dp / d6F pe = -A6 x ( gp+p _)( - — ) = , (2.9) 

-Jo P- P+ ga 

and the magnetic flux only depends on a and k and does not on p + . Thus the size of the bulk is not fixed by 
flux conservation. To discuss the modulus dynamics, we take the moduli approximation, namely assuming that 
p + — * p + (x^). Integrating over the extra dimensions, we obtain 

After redefining the modulus as 



' 8ttkM 6 4 



X6(^) = W^ JL P+, (2-11) 



we obtain the canonical form of the modulus kinetic term as 

i-^d Xe ) 2 )- (2.12) 



f.S,). / d\~r( ' 



III. THE ONE-LOOP EFFECTIVE POTENTIAL OF THE VOLUME MODULUS 

Next, we introduce a massless, minimally coupled scalar field and work in the Euclideanized space. 4 The action 
for the massless scalar field perturbations is given by 

^scalar = ^ / d? Xy/fj^Aej) . (3.1) 



4 For an earlier work on the Casimir effect in 6D spacetime, see e.g., [13 . 
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A. Scalar one-loop effective action 

The one-loop effective action for a massless minimally coupled scalar field is given by 

W 6 = ilndet(-A 6 ), (3.2) 

where Aq is 6D Laplacian. Wq needs to be regularized and renormalized. For this purpose, wc define 

,2s r oo ai , a \ 

e-* ( - Ae) ) , (3.3) 



2 J t 1 "" V 



where Tr (e ^ A|3 - 1 ) corresponds to the (integrated) heat kernel. The (integrated) zcta function is related to the heat 
kernel by a Mellin transformation: 

C(s, A 6 ) = Y^jj°° dtt^Tr ( e - ( (- A «)) = Tr((-A 6 ) s ) , (3.4) 

and after analytically continuing to s — > we obtain the renormalized one-loop effective action. The renormalized 
scalar field effective action can be written as 

W 6 , Ien = -|C'(0, A 6 ) - |C(0, A 6 )lnM 2 - (3.5) 
By integrating over the internal dimensions, the 4D effective potential is 

W 6 . rcn = f (d 4 xp 2 + )V 6 , cS = f d 4 xV 6 , cS , (3.6) 



where V e s has the dimensions (length) 4 . For brevity, from now on we shall omit the subscript "ren". 
The zeta function is given by the summation 

C(«,A 6 ) = /^£/|^L (3.7, 

where the eigenvalues are defined by 

A 6 <K = -A 2 A . (3.8) 

It is straightforward to show that 

C(0,A 6 ) = a 6 (/ = l), (3.9) 
where ag(f) is a heat kernel coefficient, defined by the asymptotic expansion of the heat kernel (l3l. [l7j: 

Tr (e-*(- A6 )) ~ Y, t {k ~ 6)/2 a k (f) , t . (3.10) 

fc>0 

B. Continuous conformal transformations 

One strategy to evaluate the one-loop effective action and the effective potential is to define a continuous conformal 
transformation (parameterized by e) 



dsl = e 



2 ^-V"dsl, W = ±ln(2p), (3.11) 
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and thus 
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where for e = 1 we have the original metric, which we shall denote as A^^ = Aq. The classical action of this scalar 
field is changed under a conformal transformation Eq. (|3.11[) 



5' 



scala 



d 6 x^gcj)A 6 



d 6 x^(A 6 + E 6 (e))^ 



(3.13) 



where 



E 6 (e) = -4(e- I) 2 g ab V a cuV b uj + 2(e - l)A 6 lncj 
2p) (l + a 2 z 2 ) 2 



(3.14) 



We perform zeta function regularization explicitly because in the unwarped frame e = we can derive exact mass 
spectrum and thus perform exact summing up of all the relevant KK modes (See Sec. IV). The correction associated 
with such a conformal transformation is commonly known as the cocycle function: 



W 6 = -ic'(0,A 6 )-^C(0,A 6 )l n/ i 2 



(3.15) 



= -ic(0,A 6 , e=0 )ln M 2 + {-ic'(0,A 6:e=0 )~^ dea e (f = d e lnQ e )}. 
The term is given by the volume integration of linear combinations of cubic order curvature invariants (l3l . Il7j : 
<«(}(/) := {^)~ 3 {J d 6 x^g (l&R'V*;* + 17R, A R A - 2R AB . C R AB ' C - AR AB . C R AC - B + 9R ABC d-eR ABCD]E 



I M 

-28RR fA ' A — 8R AB R AB; c' C + 2AR AB R A C -bc + ^RabcdR ABCD ' E ,e + -jrB 3 — —RRabR 

y 



14 

y 



AB 



14 
~~3 

-t-RabcdR^" efR l 
y 



' ' pp pABCD 2 08 n pA pBC | l "*p p pACBD lu p pA pBCDE 
—KK A BCDtt f—tiAB-ti + —i^ABi^CDn- —n A sn CDE-K 

y 



64 

y 



16 

y 



80 

—RabcdR a e^fR l 
y 



^((i/'.V + 60£ 6 £V% + 30£f£ 6; .4 + 60£| + 10i?£VV + AR AB E & , AB + 12B A E 6 , A 



J_ 

360 

+30E%R + 12E 6 B A , A + 5E 6 R 2 - 2E 6 R AB R AB + 2E 6 R ABCD R ABC d 
+ (contribution of conical branes) j . 

It is rather useful to use the effective potential, as given by Eq. (|3.6p . 

A 6 (a,K,g) - B 6 (a,n,g) ln(p 2 p + ) 



V 6t es{a,K,g,fi;p+) = 



(3.16) 



(3.17) 



where we define 



d 4 xAg(a, k, g) 



J d 4 xBQ(a,K,g) = J 



d 4~ A 6 (a,K,g) _ 
p\ 

4 „B 6 (a,K,g) 1 . . 

d x g = -C(0, A 6i£=0 ) ■ 

P+ * 



dea 6 (/ = ±ln(^)) 
2 p+ 



-C'(0,A 6)e=0 ), 



Clearly, if B§(a, n,g) > 0, then the modulus effective potential has a minimum at 

p * + _ A1 -2 e (2A 6 + B 6 )/(2B 6 ) _ 

After a redefinition of the modulus, as given by Eq. (|2.1ip , the effective potential can be rewritten as 

( 2 2 2 \ 



V^ ef f(a,K 1 g 1 p;p+) 
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xi 



The field value at the minimum is then given by 



X6,* 



8ttkM 6 4 (2A 6 +b 6 )/(2B 6 ) 
MS 



(3.18) 



(3.19) 



(3.20) 



(3.21) 
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C. Phenomenological implications after volume stabilization 

1. For the hierarchy problem 

One of the most longstanding problems in phenomenology is the hierarchy problem, namely why gravity is so weak 
in comparison to the electroweak interaction (why the Planck scale Mpi ~ 10 19 GeV is much larger than that of the 
electroweak scale Mew ~ 10 3 GeV). A way to resolve the hierarchy problem in braneworld set-up was first proposed in 
the large extra dimension model given in Ref. [2] . The basic idea is that we assume that the fundamental gravitational 
scale is not Mpi itself but the higher-dimensional one {Mq in 6D braneworld) and Mq ~ Mew- Then, the observed 
Planck scale is effectively given as a result of a dimensional reduction and in the present model 

M a ^ P±M M 4 (3 22) 

1 /2 

To get the observed value of reduced Planck scale, the size of extra dimension should be (p+) " ~ 0.1mm. Thus now 
we ask whether volume stabilization at this scale is available in the present model. 

If we assume a brane localized field, e.g., a Higgs field, 5 whose bare mass is given by m 2 on either brane at p± 
then the observed mass scales are 

m 2 + = rn 2 , m 2 _ = a 2 m 2 . (3.23) 
We now assume that m\ ~ M§ w . Thus, the mass ratio between the field and the effective Planck mass is given by 

m± _ / M 2 m 2 \ {2 A e+ B 6 )/(2B 6 ) „ / M 2 m 2 \ (g 2 a 2 \ (2Ati+Bfi)/(2Be) ( . 

M 2 ~ V Ml )\2wJ ' M 2 ~ V Ml A 2mJ ' 1 ' 

Assuming that the factor of (pm/M§) 2 takes the optimal value of 0(1) for the unification of all the fundamental 
energy scales in 6D, the effective mass ratio is characterized by 

R ( a ^ :g ) : =jL e -^ 6 +B 6 )/(2B 6 ) ^ (3 25) 

where we have used the value of given by Eq. (|3.19p . As is explained above, once the size of the internal space 

p 1 / 2 is stabilized at 0.1mm, then R has a value as ~ 10~ 32 and our main task is to explore such a possibility in the 
parameter space (a, k, g). 

2. For the cosmological constant problem 

Another important problem is the cosmological constant problem, namely why the energy density of the present 
dark energy, assuming that its origin is vacuum energy of quantum fields, pvac ~ (10 _3 eV) 4 is much smaller than Mpj, 
expected from the standard model. A characteristic property of 6D braneworld is that the tree level vacuum energy 
of the brane, i.e., brane tension, only changes the bulk deficit angle and does not affect the brane geometry Note 
that there are several criticisms for such a self-tuning mechanism mainly because flux conservation/quantization may 
induce a hidden fine-tuning of brane tension ||. Here we focus on how important the Casimir effect is at one- loop 
order of the vacuum energy, assuming that a self-tuning exists. Of course, to build more realistic models with such a 
self-tuing mechanism is an important subject, but is out of scope of this article. 

After volume stabilization, from Eqs. (|3.17p and (|3.19[) , the effective potential (Casimir energy density) takes the 
value 



V, 



r 6 V(a, «, g) = -^ 4 B 6 (a, k, g ) e -{2A e +B e )/B e (= ^ ( (3 _ 26) 



!>+=!>+., 



and hence, the realized brane vacuum energy is almost completely determined by the renormalization. The renormal- 
ization scale p is also a free parameter but in order for an additional hierarchy problem not to appear it should be 
set to a value such as \i ~ M 6 (~ Mew)- Then, B 6 oc g A and thus we get the relation depicted Eq. (jl.ip with n = 2. 
Note that Casimir energy density becomes negative as in Eq. (|3.26[) . Thus, some additional mechanisms to uplift the 
potential minimum must exist. In this article, however, we just focus on the absolute value of the energy density. 



5 It is difficult to treat any matter on a strictly codimension two brane, unless we regularize the brane. Thus, we implicitly assume that 
the brane has a small, but finite thickness rather than existing as a strictly codimension two defect. 
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IV. EVALUATION OF THE DERIVATIVE OF ZETA FUNCTION IN THE UNWARPED FRAME 

A. Mass spectrum 
From Eq. (|3.13|) . the mass spectrum in the unwarped frame is determined by 

(A 6 + £ 6 (O))0a = -A 2 0a. (4.1) 
We shall decompose the mass eigenfunction as 

/ T^E^K"^. (4.2) 
The equation of motion of equation (|4.1j) has a series of exact solutions [12, EH EH : 



1 + a 2 z 2 



ii/2k / 1 v no? /2k, 



2 



$a(^) = A/ — 5- A ;r T 2^1(1-1/ fl-aT),i/ (l-a 2 ),l- . 

w V 1 + z 2 L V! + z / V! + z / 2k v ; ' 2n y h k'1 + z 2 ' 

( z 2 \ n/2K ( 1 \ -™ q2 / 2k n „ n z 2 1 

+ B {— 2 ) (iT^j ^(1 - ^ + ^(1 - o 3 ).^ ^(1 -^.l + ^i^)]. (4-3) 

where 2^1 (&, b, c; x) is Gauss's hypergeometric function and 



1 



,2 



^^i + yi + ^ + ^i-^j. (4.4) 

Here imposing the regularities on both conical branes, we obtain the following exact mass spectrum 

9\n\ An 2 n 2 

\ 2 = k 2 + g 2 [4m(m + 1) H — — (2m + 1) (1 + a 2 ) + — — ] , m = 0,l,2,---. (4.5) 

Our method to determine the mass spectrum is essentially based on the same arguments given in [l6| (more precisely 
it is demanded that the wave functions should be Hermitian) . However, if only normalizability were imposed this 
would allow for logarithmic divergences at the poles and this may well lead to additional modes in the eigenvalue 
spectrum. But what should be stressed is that as is discussed in the case of 4D toy model [l2j the effects of these 
additional modes are not important and neglecting them is valid in the following discussions. 

B. Zeta function regularization 

The first step is the integration of the zeta function Eq. (|3.7p over the Fourier space 

(2^) 4 C( S , A 6 , £=0 ) = f d'x . ^C^^ E + 1) + — (2™ + 1)(1 + a 2 ) + ^] 2 ~ s . (4.6) 

J {S ~ W - l) ^ K K 

Hereafter, we omit the index A 6 e=0 . For a later convenience, we define 



4(1 + a 2 ) 4a 2 
a = 4, b= 1 + > , c=^-, q = -l, [3 = 1/2, (4.7) 



and using them 



-i ■ < 1+< "">h 



2a 2 2k 

q(n) = q + cn 2 - { ^- = --l(l-a 2 ) 2 -l. (4.8) 
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C. Extended binomial expansion 

Frequently, the form of the two-dimensional Epstein zeta function allows one to perform its summation in an ele gan t 
way which involves the Chowla-Selberg expansion formula or more frequently generalization of it, see [Til . [T3. Il8l|. 
Some conditions must be satisfied; the most impotant one is that the quadratic form must be positive definite and 
the constant q term should be also non-negative. But this is not here the case, and the fact that q < does not allow 
for such a beautiful analysis. Henceforth, we shall apply in what follows what we have called the extended binomial 
expansion approach, namely, we will introduce an extra summation via 



oo oo 



[a(m + fi) 2 + b(m + j3)n + cn 2 + q] 



-s+2 



m— n—1 

oo oo 



oo oo oo 



££[a(m +/ §) 2 + c s+ =££E 



r(3-s) 



m — n—1 

oo oo oo 



£ ^ (~l) j r(g + i - 2) [n(m M m 2-,- i/0 



m— n—1 j—0 

[a(m + 0) 



[a{m + P) 2 } 2 s V 



q J , 



m—0 n—1 j—0 

where the validity of the binomial expansion is defined for 

q 



a(m + (3) 2 



< 1, 



(4.9) 



(4.10) 



which is indeed satisfied for all possible values of m and n in our model. Until now, we can reduce to the following 
expression of the zeta function 



(2tt) 4 C( S ) 

/oo 
^££ 
j—0 m,n 



-2 %2g 2(2-s) r ( s + j _ 2 ) 



1 1 + a 2 



it, n 



(4.11) 



The j-summation is absolutely converging and so we can exchage the order of j-summation and (m, n)-summation. 
Then, to perform m-summation is rather straightforward. We decompose the zeta function into the contributions of 
non-axisymmetric and axisymmetric modes as 



C(*)=C(s) ^ +Co(s) 



(4.12) 



The contribution of axisymmetric modes is discussed in Appendix B and hereafter we focus on non-axisymmetric 
modes with n ^ 0. Then, we get 



(2^) 4 C(*) 



°° 00 ,o5-2s- 



^££ 

j=Q n=l 



2 j7r 2 g 2(2- s)r(s + j _ 2 ) 

~W(s) 



1 



3 t 1 1 ~\~ Q? \ 

C H (2s + 2j-4,- + ^— n). (4.13) 



The zeta function is convergent only Re(s) > 5/2 and thus we need to analytic continuation to s — ► 0. 

D. Zeta function regularization 



By subtracting the divergent terms from the asymptotic expansion of the zeta function (in the limit of s — > 0) and 
adding back the counterterms we obtain 



(4.14) 



(2^) 4 C( s )|^ o = |d 4 x^P( s ) + £A(j, s )G(j, S )j 



where 



oo oo 2 

P(s) := £ G(j, s) [£ { [!L)(1 - a 2 ) 2 + l] j C H (2s + 2j - 4, 



1 1 + a 2 

2 + 2k 



-n) - F 



(4.15) 
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and 



G(j ' s) = — w® — 9 



F(n,j,s) := 



-i-2s 2- 5+2s+2 i(l - a 2 ) 2j 



A(i,s) 



k/ (2s + 2j -5){l + a 2 ) 2s + 2 i- 5 

\w {a,j,s)(-) 6 + wi(a, j, s)(-) 4 + w 2 (a, j, s)(-) 2 + w 3 (a, j,s)\ , 

I K K K J 

2-B+2.+!y /c 2H-l( 1 + a 2 )5 - 2s -2 J(1 _ a 2ji« 

2s + 2j - 5 

j- w (a, Wi(a,j,g) w 2 (a,j,s) 

I ^ Cr(2s - 5) + - 4 Cfl(2s - 3) + - 2 Cr(2s - 1) 

+ t U3 (a > i 1 a)Cfl(2« + l)}. (4.16) 

E. Derivative of zeta functions 

1. Analytic continuation of the deriavative of £ — function 

We first make an analytic continuation of the subtracted zeta function. The derivative of the subtracted zeta 
function P'(0) is given in the Appendix A in detail. We show a typical example of this function in Fig. 1. The j and 
n-summations both show good convergency. 

2. Analytic continuation of counter terms 

The analytic continuation of the counterterms are given by 

d_ 

ds 



Y J (G(j,s)A(j,s)) 



OC 



d 

J=3 



(4.17) 

The result of analytic continuation of each j term is shown in Appendix A. We show an example showing the 
convergency of the summation in Fig. 2. 

F. On the heat kernel coefficients and cocycle functions 

As is explained in Sec. Ill, the coefficient Bq{ol, K,g) in the one-loop effective potential can be obtained by 

d 4 xB e (a,K,g) = ic(0,A 6 , e=0 ). (4.18) 

Here Bq(cy, K,g) includes both brane and bulk contributions. To decompose Be(a, K,g) into brane and bulk contribu- 
tions, it is useful to focus on the fact that a^(f = 1) = £(0, Aq^ = q) [13j. The bulk heat kernel is given by 

~ 4 <«*> ' T(z) 

'40320tt 2 J (1 + z 2 )(l + a 2 z 2 ) (1 + a 2 z 2 f ' 

where T{z) is a 12th order polynomial of z given in Appendix C. Note that now there is no way to evaluate the brane 
heat kernel a6,branc(/ = 1) because of a lack of mathematical formulation of conical heat kernel. Thus the pure bulk 
and brane contributions are evaluated by 

R , \ _ ng A f°° z T(z) 

^ 6 ,buikla, k, g) - 8064() ^ 2 ^ az ^ + ^ + ^ + ^ 2)6 , 

-B 6l brane(o!, k, g) = B 6 (a, k, g) - -B 6 ,buik(a, «, g) , (4.20) 
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FIG. 1: Plots for P'(0), Eq. (53), as a function of r are 
shown for k = 0.9. The red and blue curves correspond to 
truncation of the j-summation at j max = 100 and jmax = 
200, respectively. 



FIG. 2: Plots for (d/ds)(£\ GA)(s)| s _ , Eq. gT7) , as 
a function of a are shown for k = 0.9. The red and 
blue curves correspond to truncation of j-summation at 
jmax = 100 and j max = 200, respectively. 



respectively. The authors of Ref. [Tl| studied the contribution of conical branes to the total heat kernel and they 
conclude that these contributions are not negligible. We will see them in the next section. 
The cocycle function from the bulk part is given by 



f dea 6 {f = \\n{^)) 
Jo 2 v p+ 



= - / t**3 



de 



dz 



5(e,z) 



(1 + z 2 )(l + a 2 z 2 ) (l + a 2 z 2 ) 6 



In 



1 + a 2 z 2 



+ z 2 



(4.21) 



where Q(e,z) is also a 12th order polynomial given explicitly in Appendix C. Thus, we also obtain the coefficient 
^6,bulk( a > «j 9) as 



de- 



Kg 



dz ■ 



/n 80640tt 2 J (1 + z 2 )(l + a 2 z 2 ) {l + a 2 z 2 f 
Note that the total coefficient Aq[ol, K,g) is given by 



1 + z 2 



(4.22) 



/ 



d i xA 6 (a,n,g) = -^('(0, A 6 , 6 =o) + J d i xA 6 ^ hulk (a, n, g) + \ <l '.r. .!,;.!, 1;m ,,( n . /■;. ,/ i 
d 4 xA 6 (a, K,g) + / d 4 xA 6:hranc (a, k, g) , 



(4.23) 



where -Ag brane^j (?) is the brane part of the cocycle function and A$(a, k, g) is the total of the 1st and 2nd term 
of the middle step, which now we can evaluate. As the case of Be^ ranc (a, n,g), we have no direct way to evaluate 
^.branc^i K j ff)- In the following evaluations, we will use A 6 (a, k, g) instead of A e (a, k, g) to discuss phenomenological 
implications without the brane cocycle function. The evaluation of A e (a, k, g) is left for a future work. Note the brane 
contribution is also partially included in C'(0 ; ^-6.<i=o) an d thus we believe that Aq^ 



should not be dominant. 



V. VOLUME STABILIZATION AND PHENOMENOLOGICAL IMPLICATIONS 



A. Volume stabilization 



In Fig. 3, we plot i?6,bulk(a, K,g) (and thus the density of a§{f = 1)) given in Eq. (|4.20j) and the integrand of 
C(0, A 6 e= o) as functions of a for several choices of k. An important observation is that for smaller a the sign of the 
integrand of £(0, Ag e= o) becomes negative implying that the volume modulus is destabilized. In the case of 4D toy 
model discussed in Ref [llj . there is also negative brane contruibution but then the bulk effect still dominates and 
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the modulus is always stabilized. For any value of k, we obtain the critical value of a*(/«), below which the volume 
modulus is destabilized. In Fig. 4, we show the critical a* as a function of n. 



B. Implications for the hierarchy problem 

The resultant mass ratio of a brane localized field to the effective Planck scale is characterized by Eq. (|3.25[) . Our 
interest is the logarithmic scale of the ratio: 

**> = **> (£) - ( Jfeff + 5) 1o8 '° ' (5J) 

Note that we use A§ in Eq. (|4. 23(1 instead of A§. For the evaluation of the heat kernel coefficient B§{a, K,g), we use 
Eq. (|4.18[) . As we have seen, there is critical values of for each k, below which the modulus is destabilized. In 

Fig. 5, we show the plot of log 10 (i?(a, K = 0.3, <?)) as a function of a > a*(0.3) for g = 10~ 3 , 10" 1 , 10. A realistic 
value of log 10 (i?) ~ —32 is possible, but the warping parameter a should be close to a*(re) and thus, a fine-tuning is 
required though it is not as severe as the original fine-tuning. 



C. Implications for the cosmological constant problem 

We can also discuss the implications for the cosmological constant problem by evaluating the order of the effective 
energy density of the volume modulus. In Eq. (13.26|) . we made the optimal choice of [i ~ Mg ~ m. We show the plot 
of 



logio ( 



V* 6 ,cs(a, K,g) 



= log 



10 



B 6 (a,K,g) 



A 6 (a,K,g) \. 
B 6 {a, K,g) I 



(5.2) 



The degree of the ratio is somewhat sensitive to the value of the bulk dilaton coupling especially for smaller values. 
As is expected, the result is very similar to the case of the hierarchy problem. 

It appears to be possible to obtain an observationally acceptable value of the effective cosmological constant on the 
brane, but again a should be close to the critical value £**(«;), as for the previous case. In Fig. 6, we show the plot of 
logio (l^*6,cff(a,0.3,g)/^ 4 |) as a function of a > a* (0.3) for g = 10~ 3 , 10"\ 10. 



VI. CONCLUSION 



In this article, we have discussed the Casimir effect in a 6D warped flux compactification model based on a 6D 
supergravity and its implications for phenomenology, i.e., the hierarchy and cosmological constant problems. In its 
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tion of a. The red, green and blue curves correspond to 
the cases that g — 10~ 3 , 10 _1 , 10, respectively. Note that 
a > a»(« = 0.3). We take n max = 20 and j max = 100. 



4D effective theory of the model, a volume modulus appears and to stabilize the modulus we considered the Casimir 
effect induced by a massless, minimally coupled bulk scalar field. The effective potential of the volume modulus has 
the form of a Coleman- Weinberg potential. The stability itself can be determined by the sign of the coefficient in 
front of the logarithmic term Bg(a, k, g). 

There has been no mathematical formulation of the heat kernel coefficient for 6D conical branes as far as the 
author is aware of. However, we now have an exact mass spectrum in the unwarped frame and thanks to the fact 
that ae(f = 1) = £(0, A e= o)j by summing up all the modes with zeta function regularization, we can obtain the 
value of Bq. As a result, especially for a < a*(«), where a*(«;) is deficit angle dependent critical value, the volume 
modulus is destabilized because of a strong negative contribution from the brane quantum corrections. In contrast, for 
a > Q!*(k) the modulus is stabilized. After volume stabilization, we then discussed the implications for the hierarchy 
and cosmological constant problems and showed that it is possible to get observationally acceptable values of the 
ratio of effective mass scales and vacuum energy density. However, for each value of the deficit angle of the reference 
(+)-brane, the value of a should be tuned to be close to a*. 

As is mentioned in the text, in the framework of the present model it is not possible to obtain a positive vacuum 
energy density and this fact requires some modification of the present model, e.g., some kind of uplifting mechanism 
and/or field content contributions to the Casimir effect from the various multiplets that arise in the original super- 
gravity model. In future work, we will explore more realistic modifications of the present model. Concerning the 
background model, we need to find a 6D braneworld model where a self-tuning mechanism at tree level exists. It 
would also be intersting to consider the Casimir effect in a cosmological (time-dependent) background. We hope to 
report on these in future publications. 
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APPENDIX A: RESULTS OF ANALYTIC CONTINUATIONS 



1. Asymptotic expansion and analytically continued zeta function 



To perform analytic continuation, we derive the asymptotic expansion (See e.g., [2 
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where is the Bernoulli numbers. Following the above fact, we obtain the desired asymptotic expansion 
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Thus, we obtain the analytically continued result of the derivative of the subtracted zeta function given by Eq. (|4.15|) 
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Next we focus on the quantities appearing in the counter terms. The analytic continuation of the j = term is 
given by 



A( G ( 0)fl )A(0,-)) 



k(1 + a 2 ) 5 7r 2 <? 4 



10 



2 In 



uiq(q;, 0,0) wi(a,0, 0) u> 2 (a,0,0) 



120k 4 



ff 2 (l + a 2 )) + lo) ( 2.~>2/+ 

1 uV(a,0,0)) 1 w 2 '(a,0,0)) 2tu (a,0,0) , 



120 k 4 
2w 2 (a,0,0) w 



+ 



Ck("5) + 



12k 2 



2tui(a,0,0) , 



Ck(-3) 



where we have used 



w 3 (a,0,s)( R (2s + l) 



C R (-1) + d s (w 3 (a, 0, S )C«(2 S + l))^ Q 



3 1 _ 31 (¥- 307) f , 2) 
126(1 + a 2 ) 6 1890 (1 + a 2 ) 6 



(A5) 



(A6) 



Similarly, contributions from j = 1 and j = 2 are given by respectively 



ds 



(g(m)a(i,*)) 



K(l + a 2 ) 3 (l-a 2 ) 2 7rV 



Wo(a, 1,0) wi(a, 1,0) w 2 (a,l,0) 



120k 4 



+ 
+ 
+ 



( 21ll (.g 2 (l + a 2 )) + 3X 252k6 
( W3 ( a ,l, S )C(2 S + l))_ ) 
1 Wl '(a,l,0)) 1 w 2 '(a,l,0)) 2ioo(a,l,0) w 



120 k 4 
2w 2 (a,l,0) 



12 K 2 



+ 



-&(-5) + 



12k 2 



2wi(a,l,0) 



Ck(-3) 



where 



w 3 {a,l,s)C R (2s + l) = 



C R (-l) + d a (w3(a,l,8)Oi(28 + l)) 



-(29 + 89 a 2 + 29 a 4 ) 



(A7) 



315(-l + a 2 ) 2 (l + a 2 ) 6 

(1555 + 3358 a 2 + 1555 a 4 -24 7 (29 + 89 a 2 + 29 a 4 )) s 
3780 (-l + a 2 ) 2 (l + a 2 ) 6 



0( S 2 ). 



(A8) 



and 



K(l + a 2 )(l-a 2 ) 4 7r 2 5 4 



+ 



( 21n (^(Tf^)) + 2 )(- 252k« 
( W3 (a,2,0)C(2 S + l))^ Q ) 



iu (a,2,0) ioi(a,2,0) tu 2 (a,2,0) 



+ 



120k 4 



12k 2 



16 



1 Wl '(a,2,0)) 1 w 2 '(a,2,0)) 2w (a,2,0) ^ 2^ (a, 2,0) 

" 120 k 4 ^ ^ + ^ C fl ( 5 )+ ^ C*( ^ 



+ 2WM, 2 2 ' 0) Ca(-l) + d s (w 3 (a, 2, s )( R (2s + 1)) J , (A9) 



where 



+ 



w 3 (a,2, S )Cit(2s + 1) 
87 + 296 a 2 + 914 a 4 + 296 a 6 + 87 a 8 
630(-l + a 2 ) 4 (l + a 2 ) 6 

(-1151 + 10447 + 4 (-1663 + 8887) a 2 + 6 (-759+ 18287) a 4 + 4 (-1663 + 8887) a 6 + (-1151 + 1044 7 ) a 8 ) s 

3780 (-l + a 2 ) 4 (l + a 2 ) 6 

+ 0(s 2 ). (A10) 
Finally, for the contribution from j > 3, we obtain 
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Applying similar scheme of analytic continuation, we obtain 
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APPENDIX B: ZETA FUNCTIONS FOR AXISYMMETRIC KALUZA-KLEIN MODES 

We now derive the derivative of zeta functions for axisymmetric modes. The axisymmetric zeta function can be 
written as 

f ~ 2 5-2 8 -2j 2 2(2-s)pf • , _ 0) 1 

( 2 .) 4 co( S ) = y - — ms ) Ch(2s + 2j - 4 >2 ] - (B1) 

After taking derivative with respect to s, analytic continuation to s — > gives 

(27r)%'(0) = y d 4 x ,g 4 ( - 8( R (3) + 24^) , (B2) 



where we have used 
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and Cr( z ) = 7I ' Z ^ 2 Cr(1 ~ z )r((l — z )/2)/r(z/2). The analytic continuation of zeta function itself is 

(2^) 4 Co(0) = 0. (B4) 

APPENDIX C: FUNCTIONS RELATED TO THE HEAT KERNEL COEFFICIENT 

The function T(z) in Eq. (|4.19p is defined by 

T{z) := 1748 + 4964 z 2 + 5461 z 4 + 2035 z 6 - 1740 a 2 - 5936 z 2 a 2 - 9049 z 4 a 2 - 2882 z 6 a 2 

+ 711 z 8 a 2 + 1564 a 4 + 13944 z 2 a 4 + 32534 z 4 a 4 + 26109 z 6 a 4 + 12757 z 8 a 4 

+ 3652 z 10 a 4 - 548 a 6 - 10480 z 2 a 6 - 27054 z 4 a 6 - 30044 z 6 a 6 - 27054 z 8 a 6 

- 10480 z 10 a 6 - 548 z 12 a 6 + 3652 z 2 a 8 + 12757 z 4 a 8 + 26109 z 6 a 8 + 32534 z 8 a 8 + 13944 z 10 a 8 

+ 1564 z 12 a 8 + 711 z 4 a 10 - 2882 z 6 a 10 - 9049 z 8 a 10 - 5936 z 10 a 10 - 1740 z 12 a 10 

+ 2035 z 6 a 12 + 5461z 8 a 12 + 4964z 10 a 12 + 1748 z 12 a 12 . (CI) 

The function Q(e, z) in Eq. (|4.21ll is also defined by 

g(e, z) := 1024 + 2580 z 2 + 2808 z 4 + 1042 z 6 + 4176 z 2 a 2 + 10176 z 4 a 2 + 11964 z 6 a 2 + 4704 z 8 a 2 

- 4392 z 2 a 4 - 15984 z 4 a 4 - 26130 z 6 a 4 - 11976 z 8 a 4 - 588 z 10 a 4 + 4368 z 2 a 6 
+ 25632 z 4 a 6 +46728 z 6 a 6 + 25632 z 8 a 6 +4368 z 10 a 6 - 588 z 2 a 8 - 11976 z 4 a 8 

- 26130 z 6 a 8 - 15984 z 8 a 8 - 4392 z 10 a 8 + 4704 z 4 a 10 + 11964 z 6 a 10 + 10176 z 8 a 10 + 4176 z 10 a 10 
+ 1042 z 6 a 12 + 2808 z 8 a 12 + 2580 z 10 a 12 + 1024 z 12 a 12 + 704 e + 2368 z 2 e + 2678 z 4 e + 1014 z 6 e 

- 1728 a 2 e - 10272 z 2 a 2 e - 19726 z 4 a 2 e - 15332 z 6 a 2 e - 4150 z 8 a 2 e + 1600 a 4 e 

+ 18912 z 2 a 4 e + 50460 z 4 a 4 e + 54762 z 6 a 4 e + 26062 z 8 a 4 e + 4448 z 10 a 4 e - 576 a 6 e 

- 15456 z 2 a 6 e - 55324 z 4 a 6 e - 80888 z 6 a 6 e - 55324 z 8 a 6 e - 15456 z 10 a 6 e - 576 z 12 a 6 e 
+ 4448 z 2 a 8 e + 26062 z 4 a 8 e + 54762 z 6 a 8 e + 50460 z 8 a 8 e + 18912 z 10 a 8 e + 1600 z 12 a 8 e 

- 4150 z 4 a 10 e - 15332 z 6 a 10 e - 19726 z 8 a 10 e - 10272 z 10 a 10 e - 1728 z 12 a 10 e + 1014 z 6 a 12 e 
+ 2678 z 8 a 12 e + 2368 z 10 a 12 e + 704 z 12 a 12 e + 32 e 2 + 152 z 2 e 2 + 204 z 4 e 2 + 84 z 6 e 2 

- 48 a 2 e 2 - 464 z 2 a 2 e 2 - 1028 z 4 a 2 e 2 - 872 z 6 a 2 e 2 - 260 z 8 a 2 e 2 + 480 z 2 a 4 e 2 

+ 1688 z 4 a 4 e 2 + 1964 z 6 a 4 e 2 + 764 z 8 a 4 e 2 + 8 z 10 a 4 e 2 + 16 a 6 e 2 - 176 z 2 a 6 e 2 - 1368 z 4 a 6 e 2 - 2352 z 6 a 6 

- 1368 z 8 a 6 e 2 - 176 z 10 a 6 e 2 + 16 z 12 a 6 e 2 + 8 z 2 a 8 e 2 + 764 z 4 a 8 e 2 + 1964 z 6 a 8 e 2 + 1688 z 8 a 8 e 2 
+ 480 z 10 a 8 e 2 - 260 z 4 a 10 e 2 - 872 z 6 a 10 e 2 - 1028 z 8 a 10 e 2 - 464 z 10 a 10 e 2 - 48 z 12 a 10 e 2 

+ 84 z 6 a 12 e 2 + 204 z 8 a 12 e 2 + 152 z 10 a 12 e 2 + 32 z 12 a 12 e 2 - 12 e 3 - 132 z 2 e 3 - 222 z 4 e 3 - 102 z 6 e 3 + 36 a 2 e 

+ 608 z 2 a 2 e 3 + 1490 z 4 a 2 e 3 + 1324 z 6 a 2 e 3 + 406 z 8 a 2 e 3 - 36 a 4 e 3 - 1032 z 2 a 4 e 3 - 3544 z 4 a 4 e 3 

- 4378 z 6 a 4 e 3 - 2042 z 8 a 4 e 3 - 212 z 10 a 4 e 3 + 12 a 6 e 3 + 768 z 2 a 6 e 3 + 3912 z 4 a 6 e 3 

+ 6312 z 6 a 6 e 3 + 3912 z 8 a 6 e 3 + 768 z 10 a 6 e 3 + 12 z 12 a 6 e 3 - 212 z 2 a 8 e 3 - 2042 z 4 a 8 e 3 - 4378 z 6 a 8 e 3 

- 3544 z 8 a 8 e 3 - 1032 z 10 a 8 e 3 - 36 z 12 a 8 e 3 + 406 z 4 a 10 e 3 + 1324 z 6 a 10 e 3 + 1490 z 8 a 10 e 3 + 608 z 10 a 10 e 3 
+ 36 z 12 a 10 e 3 ~ 102 z 6 a 12 e 3 - 222 z 8 a 12 e 3 - 132 z 10 a 12 e 3 - 12 z 12 a 12 e 3 - 4 z 2 e 4 - 10 z 4 e 4 - 6 z 6 e 4 

+ 16 z 2 a 2 e 4 54 z 4 a 2 e 4 + 52 z 6 a 2 e 4 + 14 z 8 a 2 e 4 - 24 z 2 a 4 e 4 - 116 z 4 a 4 e 4 - 154 z 6 a 4 e 4 

- 66 z 8 a 4 e 4 - 4 z 10 a 4 e 4 + 16 z 2 a 6 e 4 + 124 z 4 a 6 e 4 + 216 z 6 a 6 e 4 

+ 124 z 8 a 6 e 4 + 16 z 10 a 6 e 4 - 4 z 2 a 8 e 4 - 66 z 4 a 8 e 4 - 154 z 6 a 8 e 4 - 116 z 8 a 8 e 4 - 24 z 10 a 8 e 4 + 14 z 4 a 10 e 4 

+ 52 z 6 a 10 e 4 + 54 z 8 a 10 e 4 + 16 z 10 a 10 e 4 - 6 z 6 a 12 e 4 - 10 z 8 a 12 e 4 - 4 z 10 a 12 e 4 + 3 z 4 e 5 

+ 3 z 6 e 5 - 15 z 4 a 2 e 5 - 18 z 6 a 2 e 5 - 3 z 8 a 2 e 5 + 30 z 4 a 4 e 5 + 45 z 6 a 4 e 5 + 15 z 8 a 4 e 5 - 30 z 4 a 6 e 5 

- 60 z 6 a 6 e 5 - 30 z 8 a 6 e 5 + 15 z 4 a 8 e 5 + 45 z 6 a 8 e 5 + 30 z 8 a 8 e 5 - 3 z 4 a 10 e 5 

- 18z 6 a 10 e 5 -15z 8 a 10 e 5 + 3z 6 a 12 e 5 + 3z 8 a 12 e 5 . (C 
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